arXiv:1507.03555v5 [physics.optics] 26 Feb 2016 


Necessity of Time-Reversal Symmetry Breaking for the Polar Kerr Effect in Linear 

Response 

Weejee Cho and Steven A. Kivelson 
Department of Physics, Stanford University, Stanford, California 94305, USA 
(Dated: February 29, 2016) 

We show that, measured in a backscattering geometry, the polar Kerr effect is absent if the nonlo¬ 
cal electromagnetic response function respects Onsager symmetry, characteristic of thermodynamic 
states that preserve time-reversal symmetry. A key element is an expression for the reflectivity 
tensor in terms of the retarded Green’s function. 
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The polar Kerr effect [1] refers to rotation of the 
polarization of light upon reflection due to a compo¬ 
nent of the magnetization perpendicular to the reflect¬ 
ing surface. More generally, it is interpreted as a mea¬ 
sure of the corresponding pattern of broken symmetry. 
Inevitably, a sharp onset of Kerr rotation indicates a 
symmetry-breaking phase transition. Observation of this 
phenomenon has, indeed, played a crucial role in identi¬ 
fying broken time-reversal symmetry in unconventional 
superconducting states, e.g., of Sr 2 Ru 04 [2], UPta [3], 
and URu 2 Si 2 [4]. Kerr onsets have also been found in the 
pseudogap regime of the cuprates [5-8] , although they are 
somewhat rounded and their interpretation less clear. 

In Refs. [2-8], the polar Kerr effect was measured 
in a backscattering geometry. In the linear response 
regime, the corresponding Kerr rotation originates from 
the ac Hall conductivity in the dimensions perpendicular 
to the wave vector of light [9, 10], provided that spa¬ 
tial dispersion (nonlocality) is negligible [11]. For this 
Hall conductivity to be nonzero, the medium must break 
time-reversal symmetry and mirror symmetries about all 
planes parallel to the wave vector [12]. 

However, given the extremely high sensitivity of the 
Kerr measurements (< lOOnrad), it is conceivable that 
a detectable Kerr rotation arises entirely due to spatial 
dispersion effects in a medium with zero Hall conduc¬ 
tivity. An important question is whether the symmetry 
requirements discussed in the previous paragraph carry 
over even if spatial dispersion is taken into account. It 
is easy to see that the mirror symmetries must be bro¬ 
ken in any case. The sense of Kerr rotation is reversed 
in a mirror parallel to the wave vector, so if the system 
is invariant under this mirror reflection, the Kerr angle 
must vanish. Consequences of time-reversal symmetry 
are more subtle. Kerr rotation is typically described by 
the macroscopic Maxwell equations which are not invari¬ 
ant under time reversal in the presence of dissipation, 
even if there is no broken time-reversal symmetry per se. 

Indeed, there have been proposals [13, 14] that the po¬ 
lar Kerr effect in backscattering can result from natural 
optical activity [15], a spatial dispersion effect, even if 
time-reversal symmetry is unbroken. In a medium with 
natural optical activity, the speed and damping of circu¬ 


larly polarized light depend on the handedness, and it is 
plausible that this gives rise to nonzero Kerr rotation. In 
Refs. [13, 14], the authors computed the Kerr angle in the 
long-wavelength limit, where spatial dispersion is treated 
to first order in the wave vector, and obtained nonzero 
results. Several researchers [16-19] recently adopted this 
idea to interpret Kerr signals from unconventional super¬ 
conductors. 

Yet there were earlier studies suggesting the contrary 
[20, 21]. They noted that the macroscopic Maxwell equa¬ 
tions should be consistent with the Onsager symmetry 
[22] of the electromagnetic response function. This is 
a consequence of time-reversal symmetry and thermal 
equilibrium, and holds whether the medium is dissipa¬ 
tive or not. In the long-wavelength limit, this consider¬ 
ation leads to electromagnetic boundary conditions (at 
the boundary of the reflecting medium) different from 
the ones used in the studies that found nonzero Kerr 
rotation [13, 14]. The Kerr angle computed with these 
modified boundary conditions vanishes. Moreover, there 
were arguments that did not rely on the long-wavelength 
approximation but reached the same conclusion [23, 24]. 

The significance of Onsager symmetry has been re¬ 
emphasized in a number of recent papers [25-28]. This 
has led to the retraction [29-31] of the proposals that 
various measured Kerr signals are due to optical activity 
alone. 

Still, as far as we know, the consequences of Onsager 
symmetry have not been fully clarified. One should be 
able to see how it constrains the reflectivity tensor in such 
a way that Kerr rotation is forbidden. However, existing 
results, if not relying on the long-wavelength limit, either 
bypass dealing with the reflectivity tensor [23, 26] or take 
the principle of optical reciprocity [32, 33] as the funda¬ 
mental postulate [24, 25], but this principle really should 
be derived from the Onsager symmetry of the response 
functions. The main obstacle is that without an approx¬ 
imation, there is no obvious way to express the reflec¬ 
tivity tensor in terms of the (nonlocal) response function 
of the scattering medium. In this work, we avoid this 
difficulty by expressing the reflectivity tensor in terms 
of the retarded Green’s function of the electromagnetic 
wave equation. Then, it can be shown that the symmetry 
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of the response function is inherited by the Green’s func¬ 
tion and, hence, by the reflectivity tensor. From this, we 
demonstrate, in the framework of nonlocal electrodynam¬ 
ics, that Onsager symmetry leads to optical reciprocity 
in reflection [34] and, as a special case, the absence of the 
polar Kerr effect in backscattering. 

We begin by considering the problem of light reflection 
from an arbitrary linear scattering medium. Without loss 
of generality, we assume that the medium is entirely in 
the right half-space (z > 0), as in the example illustrated 
by Fig. 1. Suppose that a right-moving plane wave of 
frequency uj and transverse wave vector k|| = {k'^,ky) 
is incident on the medium. This wave is of the form 
gzk^ r harmonic time dependence is assumed 

throughout.) Here, k'|_ = (k]|,-|-fc^) with [35] 

fsgn(u;)Y^w2/c2 -k[|2 (mi < |w|/c), 

K = I _ ( 1 ) 

(|k'| > |w|/c), 

where c is the speed of light in vacuum; £j^ can be any 
vector satisfying k'^ • = 0. One should seek a solution 

to the electromagnetic wave equation with the following 
properties: for z < 0, the solution is the sum of the inci¬ 
dent wave 6*’^+ '’ £± and a left-moving term corresponding 
to the reflected wave; for z > 0, it represents the trans¬ 
mitted wave. 

The wave equation is linear, and so is the relationship 
between the incident and reflected waves. The general 
form of a reflected (left-moving) wave linearly related to 
is 

/^e*‘^-'R.(k,„k;)-5, (z<0). (2) 

Here, k_ = (ky, —k^), and k^ is defined analogously to k'^. 
The reflectivity tensor Rtj(k||,k||) satisfies the transver- 

sality conditions k_ • Rt^(k||,k') = R(^(k||,k') • = 0. 

The backscattered component of the reflected wave is 
given by the integrand of Eq. (2) at ky = —k'. Its 
polarization can be rotated relative to that of the incident 
wave only if the antisymmetric part of R^^ (—k^, k||) does 
not vanish. That is, Kerr rotation is absent when 

R.(-k;,k;) = [R„(-k;„k;)]^. ( 3 ) 

If the medium is uniform in x and y, the reflectivity ten¬ 
sor has the form 

R4k„,k;) = (2^)25(2)(k,, (4) 

In this case, Eq. (3) reduces to a simpler condition 
R^(0) = [R.(0)]^. 

The next step is to relate light reflection to the re¬ 
tarded Green’s function in the frequency domain, de¬ 
noted by Gt^(r, r'). It is a tensor such that Gt^(r, r') ■ ^ 


z = z' z = 0 



FIG. 1. An example illustrating the relation between the 
Green’s function and the reflectivity tensor: The source 
J(r) oc 5{z — z') ^, where T z, generates a pair of counter- 
propagating plane waves. One of them is reflected from the 
medium. 


is the electric field at r due to a harmonic point source 
iw/j,oJ(r) = — r')^, where | is an arbitrary unit 

vector. In the presence of a medium, G;^ also includes 
the effect of scattering. One can decompose the elec¬ 
tric field into the incident wave G° (r, r') • ^ and the re¬ 
flected wave [Gij(r,r') — G° (r,r')] ■ where G° is the 
retarded Green’s function for free space, provided that 
both r and r' are in the half-space not containing the 
medium (z,z' < 0). This observation enables us to de¬ 
scribe light reflection using the Green’s function. How¬ 
ever, the incident wave here, generated by a point source, 
is a spherical wave. To make a connection to the reflec¬ 
tivity tensor, it is essential to construct an object that 
produces an incident plane wave. 

Such an object is, in fact, well known. Consider the 
following representation of G° [36] : 



( 2^)2 


ikbCry-r;,) 



( 5 ) 


with 



j ik'^\z-z'\ 2 

^^P(k')--.(z-zOzz. (6) 


Here, -I- and — are taken for z > z' and z < z', and 
'P(k) = I — fc“^kk. Thus, the field generated by the 
planar source distribution iujfio3{r) = e*'^i'''"(5(z — z')| is 




This corresponds to a pair of plane waves with the wave 
vectors k(^ propagating in the two half-spaces z > z' and 
z < z'. 

Now suppose that a scattering medium is placed in one 
of the half-spaces. We assume that z' < 0 and, as before, 
that the medium is in the region z > 0 (e.g., as in Fig. 
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1). Among the two plane waves, only the one in the half¬ 
space z > z' \s reflected from the medium. For z' < z < 
0, the right-hand side of Eq. (7) reduces to with 

’^(k+) ■ k- The resulting reflected 
wave can be either obtained from Eq. (2) or expressed 
in terms of the Green’s function. Thus, for z,z' < 0, 


le 


-ik'z' 


2k' 


(27r)2 


AVi- r 




= I d2r;e*‘^:-ri[g^(r,r')-G°(r,r')] -I 


( 8 ) 


Here, -I- and — are taken for z > z' and z < z', 
and the transversality conditions read k± • T*(k||,k') = 

T*(k||,k||) • k^ =0. Eqs. (12) and (14) lead to the 
reciprocity relation for transmission [34]: 

lTi(k,„k;) = l[T^(-k;„-k„)]^. (15) 

We now derive Eq. (12). Consider general linear con¬ 
stitutive relations for time-harmonic fields in nonlocal 
electrodynamics: 


and Rt„(k||,k') ■ 'P(k'^) = R,^(k||,k') due to transver¬ 
sality. As ^ is arbitrary, Eq. (8) without the contrac¬ 
tion with 4 holds as a tensor identity. Then, the Eourier 
transform with respect to ry allows us to express the re¬ 
flectivity tensor in terms of the Green’s function as 


2fcC 


e-*('=^"+'=^"')R.(ky,k;) 


(9) 


= 0c.;(z,z';k||,ky) - 0°(z,z';k||,ky) {z,z' < 0), 

where 

g.(z,z';k||,k;)= / d\dA\e-d^r’^^^AAGUr,r') 

( 10 ) 

and 


D=i<,E, H = ^oiB, (16) 

where the permittivity operator is defined by 

l(?<,E)(r)=E(r)+ [ dA ■ E{r'). (17) 

Co J 

We have adopted the Landau-Lifshitz approach (§103 of 
[37]), in which the effect of the medium is solely con¬ 
tained in 'euj, i.e., the permeability is simply taken to be 
the constant /io- The kernel Xui is the electromagnetic 
response function (nonlocal susceptibility tensor). Erom 
the linear response theory, it can be shown that time- 
reversal symmetry and thermal equilibrium lead to the 
following Onsager symmetry relation [22]: 

X<.,(r,r') = [Xa;(i’»]^- (18) 


0° (z, z'; ky, kj,) = ( 2 ^) 2 ^( 2 )(k^^ _ k;) ^0 

As will be shown later, the Onsager symmetry of the 
electromagnetic response function implies 

e,^(z,z';ky,ky) = [a^(z',z;-k',-ky)]^. (12) 

Eqs. (6) and (11) show that also satisfies the analo¬ 
gous relation. (This can be viewed as a special case of the 
above equation.) From Eqs. (9) and (12), the reciprocity 
relation for light reflection follows: 

^R.(ky,k;) = l[R^(-k;,-ky)]^, (13) 

which was previously proved under the assumption that 
spatial dispersion effects are negligible [34]. Eq. (13) 
reduces to Eq. (3) for backscattering (ky = — k') and, 
hence, implies the absence of Kerr rotation. 

As a side note, a similar construction can also be ap¬ 
plied to light transmission. If the scattering medium has 
a finite extent in the z direction, one can define, in a man¬ 
ner analogous to the reflectivity tensor, the transmissiv¬ 
ity tensors for right- and left-moving waves (denoted as 
T*). It is straightforward to verify that if z and z' lie 
outside of and on opposite sides of the medium. 


2k', 



(k|hk]| 


) = 9u{z,z']'k\^,'k\^. (14) 


The dynamics of the fields is given by the macroscopic 
Maxwell equations: 

V X E = iujB, V • B = 0, 

19) 

V-D=p, VxH = J- iwD. ^ ^ 

Eqs. (16) and (19) lead to an integro-differential equation 
for E, the electromagnetic wave equation: 

TlcijE = (w^o) -f iwe^jE = J, (20) 

where At is the transverse Laplacian defined by A^E = 
—V X (V X E). Recall that we assume the scatter¬ 
ing medium is entirely in the region z > 0; that is, 
^^(rjrO = 0 if z < 0 or z' < 0 [38]. Hence, 

{C^E){r) = ^(At + ‘4)nr) {z<0). (21) 

V c J 

The Green’s function is the kernel associated with the 
operator (iuj: 

-^{C-J3){r) = [ d3r'G.(r,r') ■ J(r'). (22) 

This definition shows that the symmetry of Ga;(r,r') is 
inherited from that of In particular, is a complex 
symmetric operator, i.e., 

(Ei,£tdE2) = {£tjEi,E2), 


(23) 
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for El and E 2 belonging to a suitable class of vector¬ 
valued functions. Here, we define 



Eq. (23) follows from Eq. (18) and (Ei,AtE 2 ) = 
(AtEi, E 2 ) (integration by parts without a surface term). 
Defining Ji = £(jEi and J 2 = £wE 2 , we see that is 
also complex symmetric: 

(£-iJi,J2) = (Ji,£:1J2). (25) 

Together with Eq. (22), this implies the symmetry rela¬ 
tion 

G<^(r,r') = [G<^(r',r)]’^, (26) 

and Eq. (12) follows by virtue of Eq. (10). 

In the above analysis, one should ensure that G(^ or, 

—1 ^ 

equivalently, is uniquely determined. Notice that G^ 
is the retarded Green’s function. In other words, C~j3 is 
an outgoing wave (with the effect of scattering included), 
meaning that no radiation originates from infinity. [Only 
then can the right-hand side of Eq. (8) be identified as a 
pure reflected wave.] Intuitively, we expect that restrict- 
ing to be the retarded Green’s function guarantees 
its uniqueness; a given current distribution should result 
in a unique outgoing electromagnetic field distribution. 

It is possible to make the uniqueness argument more 
formal. We sketch the arguments here; a more thorough 
mathematical treatment is given in the Supplemental 
Material. As is often done to ensure causality, we add an 
arbitrarily small positive imaginary part to the frequency 
(w —>■ u}^ = Lo -\- irf) [39]. The frequency shift irj makes 
the entire space slightly dissipative. This fact is readily 
seen for a homogeneous, lossless, propagating medium. 
Gonsider a plane wave proportional to and a disper¬ 
sion relation w = /(k). Upon shifting uj by irj, k should 
change by an amount dk satisfying dk ■ {df/dk) = ir] to 
maintain the dispersion relation. The additional factor 
giik r exponentially decays along the direction of propa¬ 
gation (parallel to df/dk). 

In the presence of the dissipation introduced by ir], a 
wave E(r) originating from infinitely far away must have 
a divergent amplitude there. Otherwise it vanishes at any 
finite r, i.e., after having traveled and been dissipated 
over an infinite distance. For the same reason, we expect 
that a wave is exponentially suppressed at infinity if and 
only if it is outgoing. Two such waves Ei and E 2 satisfy 
(Ei,AtE2) = (AtEi,E2) (i .e., no surface term arises) 
and, hence, Eq. (23). Moreover, one can select outgoing 
waves by demanding E to be square integrable. This 
requirement, in particular, excludes nontrivial solutions 
of the homogeneous equation £i^^E = 0 because they 
correspond to waves generated at infinity. It theiMollows 
that the solution of the inhomogeneous equation E = 
J is unique. 


In summary, we have shown, in the framework of non¬ 
local electrodynamics, that the Onsager symmetry of the 
electromagnetic response function implies the absence of 
Kerr rotation in backscattering and, more generally, the 
principle of optical reciprocity. An important observation 
is that the symmetry property of the response function 
is inherited by the Green’s function and then by the re¬ 
flectivity tensor. 
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Supplemental Material 

In the main body of the paper, a symmetry relation for 
the Green’s function [Eq. (12)] played a crucial role in 
proving the absence of the polar Kerr effect. To establish 
this relation, we argued, first, that a unique outgoing 
wave solution exists for a given source distribution and, 
second, that (Ei, AtE2) = (AtEi,E2) for two outgoing 
waves El and E 2 . [Recall that (u, v) = J d^r u(r) • v(r).] 
We formally prove these statements here. As before, the 
frequency is assumed to have an arbitrarily small positive 
imaginary part, i.e., is of the form = uj + ir]. 

To proceed further, a precise definition of an outgoing 
wave must be given. In the main text, we have examined 
wave propagation in a homogeneous, lossless medium. 
Based on this, it was argued that shifting the frequency 
by ir] leads to dissipation, enabling us to select outgoing 
waves by considering only square-integrable field configu¬ 
rations. More generally, when the medium may be inho¬ 
mogeneous, it is unclear how to explicitly describe wave 
propagation. To account for such cases, we take square 
integrability as the defining property of an outgoing wave. 
That is, given that the entire space is dissipative, we 
restrict to field configurations belonging to the Hilbert 
space 

£2 = {E : llEj] < 00 } (11E|] = ^(E^). (A.l) 

In this setting, it is important to ensure that the fre¬ 
quency shift UJ ^ UJ+ = UJ + ir] actually introduces dissi¬ 
pation. We have provided evidence for this but only in 
the case of homogeneous media. In general, one should 
demonstrate that the dissipated power [37] given by 

P^JE] = Im((E*,w,?^,E) + (H:^,u;,^oH„,)) (A.2) 

is positive, where x E. The positiv¬ 

ity condition to be verified is that 

P.,[E]>AllEf 


(A.3) 
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for some A > 0. This is different from [E] >0, which 
does not rule out the possibility that [E] vanishes for 
an “almost” square-integrable field configuration, e.g., a 
plane wave. (In such a case, Ptj^[E] can be arbitrarily 
close to zero for a fixed ||E||.) Notice that 

Pu, [E] = P^ [E] + 2r; [E], (A.4) 


where 


C/„[E] = -Re 


E*,d^E) + (H:,AioH4 


. (A.5) 


The second law of thermodynamics requires that P^ [E] > 
0 [40], where the equality holds in the absence of dissi¬ 
pation (e.g., if the medium is lossless or E is localized 
outside the medium). When Pd[E]/||E|p —0, C/cd[E] 
can be interpreted as the field energy (see Section 2.3 of 
[41]), which must be positive, i.e., C/c^[E] > kIJE]]^ for 
some K > 0. From the behaviors of Pi.d[E] and C4;[E], we 
see that Eq. (A.3) with A = 0 (? 7 ) holds. 

We now prove that the solution of the equation 

= (w+^o)”^«AtE -I- iuj+e^^'E = J (A. 6 ) 

is unique if it exists. Assuming (as we will show later to 
be the case) that there is no surface term, integration by 
parts gives (E*, — AjE) = (V x E*, V x E) and hence 

-Re(E*,£<,^E) = P^^m > AlJEJl^, (A.7) 

where the inequality is simply Eq. (A.3). Therefore, 
£ij^E = 0 implies E = 0. Suppose that there exist two 
solutions El and E 2 for a given source distribution J, 
i.e., J = = £(^^E 2 . Then we must have Ei = E 2 

because — E 2 ) = 0 . 

It remains to show that outgoing wave solutions exist 
and that one can integrate by parts without surface terms 
in Eq. (23) of the main text and in Eq. (A.7). More 
precisely, we aim to establish the following: first, that 
the equation P^j^E = J has solutions for all J in 
second, that each solution belongs to a subset of 

LP' for which the integration by parts formulas we have 
used are valid. Both statements refer to the properties of 
as a linear operator on L^. Because is a linear 
combination of and A;, we begin by considering them 
separately. 

Concerning the permittivity operator , we make an 
additional assumption for further progress: we require 
that it is a bounded operator on i.e.. 


11?.,E11 < 7llE]l (A. 8 ) 


for some 7 > 0. This condition can be derived from a 
reasonable assumption on the electromagnetic response 
function; it is sufficient if there exists a positive function 
M(r,r') satisfying 

|5ta;+(r,r') • sj < M(r,r')|sl (A.9) 


for all r, r', s in and 


J (fir' M{r, r') < (3, 



r M{r, r') < /?' 


(A.IO) 


for some /?,/?' > 0. Eqs. (A.9) and (A.IO) imply 
Ijeg^e^^E — Ejj < (/ 3 / 3 ')^A||e|| due to Schur test (The¬ 
orem 5.2 of [42]). Then, the triangle inequality, i.e., 
jju -I- vjj < |]ujj -I- [jvj] for all u, v G leads to Eq. 
(A. 8 ) with 7 = + 1]. 

The transverse Laplacian At [43] is an unbounded op¬ 
erator. The domain of At is 'D{At) = {E G : AfE G 
L^}, a dense subset of L^. This is in marked contrast 
to , whose domain defined analogously coincides with 
the entire Lfi because of Eq. (A. 8 ). A notable prop¬ 
erty of At is that it is a closed operator, i.e., if two 
sequences {E„ G X>(At)} and {AtE„} converge respec¬ 
tively to E and F both in Lfi, then E belongs to 'D(At) 
and F = AtE. [Due to Eq. (A. 8 ), {ecd+E„} converges 
to e^j^E if {E„} converges to E, so is automatically 
a closed operator.] Furthermore, At is self-adjoint, i.e., 
ViAl) = V{At) and A\ = At [44]. 

Next, consider = AAt+ where iA~^ = w+^o 
and B = ioj+, acting on V{Cuj^) = {E G : /la;^E G 
L^}. The triangle inequality and Eq. (A. 8 ) give 

ll£<.,Ell<lAl|lAtEll+7lRlllEll, 

lAlllAtEll<ll£^,Ell+7|i?|||Ell, 

which, respectively, imply that D(At) C I?(£^^) and 
that ViCujfi) C V{At). Thus, 'D{Cujfi) = D(At), and 
using the definition of a closed operator, it is straightfor¬ 
ward to check that is closed. 

The adjoint of is defined as follows: E is an element 

ofD(£iJ hand only if (E*,£^^E) = (J*, E) for a J in 

and for all E in X>(£i^^), in which case we set J = iC][,^E. 
From the fact that all bounded operators have bounded 
adjoints [44] and that At is self-adjoint, one can show 
that D(£t J = V{At) = 

We now prove that the equation J = has a so¬ 

lution for each J in That is, the range of (or, 
equivalently, the domain of C~}Z) defined by 

77(£^J = {J G : J = £^,E, E G D(£^J}, (A.12) 

coincides with Lfi. To see this, we first show that there 
is no nonzero element of orthogonal to TZ{Cuj+)- If 
such an element E^ existed, (Ej^,£tj^E) = 0 for all E 
in 'D{Ci^fi). It would then follow that Ej^ belongs to 
= 'Pifi'uj+) and that £f,^Ej_ = 0. Consequently, 
(El,A^^^Ej_) = (E;^^Zi^Ej^)* = 0, contradicting Eq. 
(A.7). Therefore, TZ{Cu>+) has a trivial orthogonal com¬ 
plement and hence is dense in L^. Furthermore, Eq. 
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(A.7) and Schwarz’s inequality lead to 
||£<,,E|| ||E|| > |(E*,£^,E)| > |Re(E*,£„,E) 

>A||Ef, 

whence 

||E = £-iJ|| <A-i||J||. (A.14) 

Eq. (A.14) implies that given a sequence {J„ € 
converging in {E„ = is also con¬ 

vergent. Let Joo and Eqo be the limits of these sequences. 
Because is a closed operator, Eoo belongs to X>(£cd+), 
and Joo = £i.j+Eoo- Therefore, Joo is an element of 
T^{^uj+), i-e., 7^(£w+) is a closed subset of As it is 
dense and closed, TZ{£ui+) = L^- 

Knowing = T*(At), we can justify 

why it was possible to integrate by parts without surface 
terms in Eq. (23) of the main text and in Eq. (A.7). It 
suffices to prove that 

(El, AtEa) = -(V X El, V X E 2 ) = (AtEi,E 2 ) (A.15) 

for all El, E 2 £ X>(A(). According to Parseval-Plancherel 
theorem, the Fourier transform defined by 

(^E)(k) = J d3rE(r)e-'^- (A.16) 

is a unitary map of onto itself. With V = E — A^E, 
we have 

||V X E|| = ||k X J-EII = 11(1 + k'^y^k X J-VII 

< 11(1 + ky-^kJ^V\\ < 2-^\\TV\\ (A.17) 

= 2-i||V|| <2-i(||E|| + ||AiE||). 

This inequality shows that if both E and A^E belong to 
L^, then so does V x E. Hence, all six quantities entering 
(•, •) in Eq. (A.15) are elements of L^. It also follows from 
Parseval-Plancherel theorem that (u,v) = ((J^u*)*, J^v) 
for all u, V S Then, it is straightforward to verify 
Eq. (A.15) by considering the Fourier transforms of the 
quantities in (•,•). 

Thus far, we have formally proved certain properties 
of outgoing waves leading to the desired symmetry of 
the Green’s function [Eq. (12) of the main text]. To be 
completely rigorous, however, one should ensure the ex¬ 
istence of the limit 77 —>■ 0+, called the limiting absorption 
principle [45, 46]. Extrapolation of our result does not 
work because, e.g., A“^ = 0{ri~y in Eq. (A.14) diverges 
in this limit. To the best of our knowledge, there is no 
result sufficiently broad in scope to establish the limit¬ 
ing absorption principle for the nonlocal electrodynamics 
problem we have considered. 
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